The RedHot analysis was completed to compare the health of western redcedar across urban heat islands. Redcedar health data were collected by community scientists in the Western Redcedar Dieback Map project on iNaturalist. Urban heat data were commissioned by three cities in the northwest using methods described in Voelkel and Shandas (2017).
Observations of western redcedar were downloaded from iNaturalist and urban heat data were downloaded from open data portals or provided by contacts in the City of Tacoma, King County (Washington) and Portland. Trees in WA were also evaluated based on EHD Ranks. HOLC data were also investigated for each city.
Note temperature data is different for each dataset and may have been collected slightly differently. Temperature data will need to be standardized (difference from mean) for each dataset, then temperatures can be compared region wide.
Environmental Health Disparity Rank Data (Full EHD Rank data (v2)) was downloaded on 6.17.23 from https://geo.wa.gov/datasets/WADOH::full-environmental-health-disparities-version-2-extract/explore
Metadata is available in this technical report: https://doh.wa.gov/sites/default/files/2022-07/311-011-EHD-Map-Tech-Report_0.pdf?uid=634dcf4aec2b5
City tree data were ‘joined by attribute’ separately so we have 3 different tree datasets to work with or merge.
Also, given the UHI data were extracted with shapefiles, the column names are limited to 10characters. Therefore, exported UHI data only include iNat ID numbers and UHI data. These were then re-merged (see below) with the iNat data to get the remaining columns with proper names.
Filter data to only include necessary columns
[1] “id”
[23] “latitude”
[24] “longitude” [35] “place_town_name”
[36] “place_county_name”
Export data for QGIS data joins
Tacoma - 357 Trees Portland - 465 Trees King County - 516 Trees
Export final .shp files as .csv
Re-import Data after following QGIS Methods
Note we needed to convert Tacoma temps to F to match king county
Some of the iNat project questions changed since it was created so some we need to adjust the answers to be more consistent throughout the project.
## Warning: Removed 36 rows containing missing values (`geom_point()`).
Should we remove these four outliers with less than -8 dist from mean af?
## Warning in geom_histogram(stat = "count"): Ignoring unknown parameters:
## `binwidth`, `bins`, and `pad`
## Warning in geom_histogram(stat = "count"): Ignoring unknown parameters:
## `binwidth`, `bins`, and `pad`
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.
The data is likely zero inflated
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.
There are some data that need to be cleaned, healthy trees should not have dieback percent of 100.
binomial distributions
Ordinal regression, maybe poisson distribution?
Some cleaning is needed, some healthy observations have dieback and some dead trees have 0% dieback.
## `geom_smooth()` using formula = 'y ~ x'
## `geom_smooth()` using formula = 'y ~ x'
Note we need to remove some outliers
## `geom_smooth()` using formula = 'y ~ x'
## `geom_smooth()` using formula = 'y ~ x'
## `geom_smooth()` using formula = 'y ~ x'
## `geom_smooth()` using formula = 'y ~ x'
## `geom_smooth()` using formula = 'y ~ x'
## `geom_smooth()` using formula = 'y ~ x'
Binary Response (2 categories)
## # A tibble: 2 × 2
## # Groups: binary.tree.canopy.symptoms [2]
## binary.tree.canopy.symptoms n
## <fct> <int>
## 1 Healthy 815
## 2 Unhealthy 424
Filtered response (5 categories)
## # A tibble: 4 × 2
## # Groups: reclassified.tree.canopy.symptoms [4]
## reclassified.tree.canopy.symptoms n
## <fct> <int>
## 1 Healthy 815
## 2 Thinning Canopy 165
## 3 Dead Top 111
## 4 Other 148
## # A tibble: 146 × 2
## # Groups: user_login [146]
## user_login n
## <chr> <int>
## 1 abbielisabeth 4
## 2 abbigail_white 5
## 3 abe 1
## 4 adrienne_stclair 2
## 5 akreiner 6
## 6 alexis_mushroom 57
## 7 alginger 4
## 8 amy_boucher 1
## 9 and_allies 3
## 10 angela_mabel 50
## # ℹ 136 more rows
## n.participants n.obs
## 1 146 1239
## n.participants n.obs
## 1 96 785
## # A tibble: 3 × 1
## Area
## <chr>
## 1 King County
## 2 Portland
## 3 Tacoma